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which  the  incoming  waves  will  encounter  as  they  enter  the  joint.  Finally,  the 
frequency  response  function  of  each  Fourier  transformed  state  variable  of  the 
transmitted  waves  is  presented  analytically . .The  results  show  that  up  to  three 
independent  waves  may  be  generated  in  the  lattice  assembly  by  externally  applied 
sinusiodal  loads:  one,  a  nondispersive  propagating  longitudinal  wave;  two,  a 
dispersive  propagating  mode  1  flexural  wave  and;  'rhree,  a  mode  2  flexural  wave. 

The  mode  2  flexural  wave  exhibits  a  characteristic  cutoff  frequency.  For  frequencies 
less  than  the  cutoff,  the  mode  2  flexural  wave  is  a  nonpropagating  spatially 
decaying,  evanescent,  wave.  For  frequencies  greater  than  the  cutoff,  the  mode  2 
flexural  wave  is  a  dispersive  propagating  wave.  From  the  scattering  matrix,  the 
"L"  lattice  junction  exhibits  dispersive  wave  scattering  and  wave  mode  conversion. 
These  features  are  illustrated  graphically  through  magnitude  and  phase  plots  of 
the  reflection  and  transmission  coefficients. 
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Introduction 

Many  design  proposals  for  large  orbital  space  systems  incorporate  lattice  type  structures. 
Modeled  as  a  network  of  joints  and  slender  beams,  such  structures  may  exhibit  low 
vibrational  damping  and  low  global  stiffness  compared  with  conventional  lattice  structures 
used  on  earth.  In  order  to  design  and  control  such  structures  effectively,  comprehensive 
dynamic  analyses  may  be  required.  The  propagation  of  structural  disturbances  through  large 
lattice  systems  is  of  practical  interest  for  design  and  nondestructive  evaluation.  Many  wave 
propagation  analyses  of  large  lattice  type  space  structures  must  consider  the  reflection  and 
transmission  of  waves  through  interconnecting  junctions.  This  report  presents  an  analysis, 
using  wave-mode  coordinates,  of  wave  propagation  through  a  two-dimensional  "L"  junction 
with  mass  and  mass  rotary  inertia. 

A  wave  propagation  analysis,  using  wave-mode  coordinates,  of  a  two-dimensional  "L" 
lattice  junction  is  presented.  The  configuration  consists  of  two  perpendicular  semi-infinite 
members  which  are  joined  by  a  rigid  joint  as  shown  in  Fig.  1 .  Each  member  is  modeled  as  a 
combined  longitudinal  rod  and  Timoshenko  beam.  The  corresponding  geometric  and  material 
properties  of  each  member  are  identical.  The  joint  is  modeled  as  a  mass  (m)  of  negligible 
geometric  extent  with  mass  rotary  inei  '  vi  J.  The  junction  is  two-dimensional  because  the 
members  attached  to  the  joint  exhibit  bou. .  dal  and  flexural  elastic  properties  such  that  the 
junction  may  translate  and  rotate  only  in  one  plane.  Points  1  and  4,  shown  in  Fig.  1,  are 
placed  equidistant  from  the  geometric  center  of  the  rigid  joint.  The  system  is  subjected  to  an 
externally  applied  sinusoidal  axial  force  TN,  shear  force  Tv  and  bending  moment  at  point  1 
(refer  to  Fig.  2a),  where  the  overbar  denotes  the  Fourier  transform  of  the  time  domain  applied 
load. 

This  report  is  presented  in  four  major  sections.  In  the  first  section,  the  concept  of  the 
joint  coupling  matrix  is  applied  to  point  1  in  order  to  determine  the  input  wave-mode  vector 
which  represents  the  waves  generated  by  the  applied  sinusoidal  loads.  The  second  section 
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introduces  the  joint  coupling  relationship  for  the  two-dimensional  rigid  joint  with  mass  and 
rotary  inertia.  The  scattering  matrix  of  the  junction  is  developed  and  presented  both 
analytically  and  graphically.  In  the  third  section,  the  results  of  the  previous  two  sections  are 
used  to  establish  the  frequency  response  function  of  each  Fourier  transformed  state  variable 
at  point  4.  Finally,  a  summary  of  the  conclusions  and  recommendations  is  presented. 

This  report  has  three  objectives:  one;  to  illustrate  the  use  of  wave-mode  coordinates  to 
study  the  propagation  of  waves  through  an  "L"  lattice  assembly  as  shown  in  Fig.  1;  two,  to 
examine  the  reflection  and  transmission  coefficients  of  longitudinal  and  flexural  waves  in  an 
"L"  lattice  assembly  at  various  frequencies  and;  three,  to  illustrate  the  phenomenon  of 
wave-mode  conversion;  that  is,  a  propagating  mode  1  flexural  wave  may  reflect  and/or 
transmit  as  an  evanescent  mode  2  flexural  wave  and  vice  versa. 
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Waves  Induced  by  Sinusoidal  Loading 
in  Continuous  Lattice  Member 

In  this  section,  the  concept  of  the  joint  coupling  matrix  is  applied  to  point  1  of  Fig.  1. 
The  wave-mode  vector  which  contains  the  waves  generated  by  the  applied  loads  at  point  1  is 
determined.  The  notation  and  definitions  presented  throughout  this  report  will  follow  those 
given  in  [1],  For  convenience,  a  brief  presentation  of  the  concepts  and  nomenclature  used  in 
wave-mode  coordinate  analysis  for  wave  propagation  studies  is  given  in  Appendix  A. 

Joint  Coupling  Matrix 

A  rigid  massless  joint  with  no  spatial  extent  is  simulated  at  point  1  as  shown  in  Fig.  2a. 
Points  L  and  R  are  located  just  to  the  left  and  just  to  the  right  of  point  1 ,  respectively. 

The  explicit  form  of  the  joint  coupling  matrix  for  arbitrary  two-dimensional  rigid  joints 
is  given  in  [2].  It  is  important  to  note  that  the  positive  orientation  of  the  internal  shear  force 
of  the  lattice  member  used  in  the  development  of  the  generalized  joint  coupling  matrix  in  [2] 
is  opposite  to  the  convention  used  in  the  derivation  of  the  transfer  matrix  of  the  lattice 
member  [3],  It  is  necessary  thajthe  two  systems  be  consistent.  Therefore,  the  polarity  of  the 
shear  force  components  of  the  joint  coupling  matrix  given  by  [2]  must  be  modified 
accordingly.  From  [2],  the  joint  coupling  relationship  for  the  rigid  joint  with  zero  mass  and 
zero  spatial  extent  shown  in  Fig.  2a  is 
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are  the  Fourier  transformed  state  vectors  at  points  L  and  R,  respectively,  and  the  superscript T 
denotes  the  transpose.  The  two  elements  of  the  joint  coupling  matrix  B  which  are  enclosed  in 
parentheses  are  the  modified  terms  which  ensure  that  the  sign  conventions  for  the  joint 
coupling  matrix  and  the  transfer  matrix  of  the  lattice  member  are  consistent. 

The  joint  coupling  relationship  contains  the  force  dynamic  and  geometric  compatibility 
requirements  at  the  joint.  Substituting  eqns.  (2)  through  (4)  into  eqn.  (1)  and  writing  out  each 
scalar  equation  give 
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Eqns.  (5)  through  (7)  represent  the  compatibility  requirements  at  the  joint  and  eqns.  (8) 
through  (10)  represent  the  force  dynamic  requirements  at  the  joint. 

Scattering  and  Wave -Mode  Generation  Matrix 

The  input-output  relationship  of  a  joint  may  be  stated  as  [2] 

wow  =  £(«& + £  (0))Lx,  0 1 ) 
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where  the  matrix  £(co)  is  called  the  scattering  matrix  of  the  joint,  and  the  matrix  £(co)  is 
called  the  wave-mode  generation  matrix  of  the  joint.  The  process  of  transforming  eqn.  (1) 
into  a  scattering  relationship  of  the  form  of  eqn.  (1 1)  is  presented  below. 

The  state  vectors  at  points  L  and  R  in  Fig.  2a  are  related  to  the  wave-mode  coordinates 
at  points  L  and  R,  respectively,  by 


where  [1,4,5] 

1 
0 


I(G»  = 


0 


0 


0 


~J ' 


AE\ 


-]■ 


-El 


0 
1 

IX 2 

(A|-q) 

I2 
.  E/p4 
'  1% 


J  /3 


0 


El 


Ll  ~  Y(<o)wE 
Lr  =  Y((i))wR 

0 

1 

fll+q) 

IX 3 

(A|+q) 
i2 
£/(34 
'  1% 

0 


(12) 

(13) 


0 

1 


flj+q) 

IX 3 


El 


flj+q) 

/2 

£/[34 


0 


0 

1 


-El 


(X|-q) 

IX, 

fll-q) 

/2 

.  EI  p4 
'  /% 


7  3 


0 


1  ^ 

0 

0 

0 

0 

4£X, 


(14) 


is  the  wave-mode  matrix  of  the  lattice  member,  and 
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wLt  =  { wL *  wL+2  vvTj  wLj  wLj  vriL,}  (15) 

wRT  =  -Tw/^T  w/?2  W??3  wR  3  w/?2  w/?-}  (16) 

are  the  wave-mode  vectors  at  points  L  and  R,  respectively.  A  summary  of  the  variables 


which  constitute  eqn.  (14)  is  given  in  Appendix  A. 
Substituting  eqns.  (12)  and  (13)  into  eqn.  (1)  gives 
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Substituting  eqn.  (14)  into  eqn.  (17)  and  multiplying  out  give 
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where  Yt.  represents  the  element  located  at  row  i  and  column  j  of  the  wave-mode  matrix  F(oo) 
given  by  eqn.  (14). 

From  Fig.  2b,  two  wave-mode  vectors  which  represent  waves  leaving  and  entering  the 
massless  joint  may  be  written  as 
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The  wave-mode  coordinate  wE*  and  wlff  represents  a  type  i  wave  (i=l,2,3)  which,  given  a 
plus  (+)  sign,  travels  along  increasing  x  (toward  the  rigid  joint  with  mass  and  rotary  inertia) 
or,  given  a  minus  (-)  sign,  travels  along  decreasing  x  (away  from  the  rigid  joint  with  mass 
and  rotary  inertia),  at  points  L  and  R,  respectively. 

Rearranging  the  order  of  the  scalar  equations  of  eqn.  (18)  gives 
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The  inverse  of  the  matrix  is  determined  using  the  computer  program  called 
MACSYMA  [6],  MACSYMA  is  a  symbolic  math  processor  that  can  perform  symbolic 
matrix  inversions  and  various  matrix  operations.  It  is  available  on  M.I.T.’s  Honeywell  DPS 


8/70M  computer  running  under  the  MULTICS  operating  system. 
Rearranging  eqn.  (22)  and  premultiplying  both  sides  by  D^m  give 
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is  the  scattering  matrix  of  the  massless  joint  at  point  1,  and 
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is  the  wave-mode  generation  matrix  of  the  massless  joint  at  point  1. 

The  upper-left  and  lower-right  3x3  matrices  of  S( co)  given  by  eqn.  (26)  contain  the 
reflection  coefficients  for  the  waves  which  strike  the  massless  joint  from  the  left  and  from  the 
right,  respectively.  The  upper-right  and  lower-left  3x3  matrices  of  5(w)  given  by  eqn.  (26) 
contain  the  transmission  coefficients  for  the  waves  which  strike  the  massless  joint  from  the 
right  and  from  the  left,  respectively.  Since  the  scattering  matrix  S(co)  given  by  eqn.  (26)  is 
independent  of  frequency,  the  system  shown  in  Fig.  2a  exhibits  nondispersive  wave 
scattering  at  point  1. 
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Consider  the  case  where  there  are  no  applied  loads  at  point  1.  Substituting  eqn.  (26)  into 


eqn.  (25)  and  multiplying  out  give 
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Eqn.  (28)  states  that  longitudinal  and  flexural  waves  may  pass  from  point  L  to  point  R 
(or  point  R  to  point  L),  through  point  1,  unaltered.  Since  each  wave  passing  through  point  1 
does  so  unaltered,  without  reflecting  or  transmitting  waves  of  a  different  type,  wave-mode 
conversion  does  not  occur  at  point  1 . 

Wave-Mode  Input  Vector 

Substituting  eqn.  (27)  into  eqn.  (25)  and  considering  only  waves  generated  at  point  1 

give 
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.here  Git  represents  the  element  located  at  row  i  and  column  j  of  the  wave-mode  generation 
matrix  £r(co)  given  by  eqn.  (27).  Eqn.  (29)  contains  the  longitudinal  and  flexural  waves 
generated  by  the  sinusoidal  axial  force  TN,  shear  force  Fv  and  bending  moment  TM  at  point  1 . 
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Joint  Coupling  and  Scattering 
Matrix  of  "L"  Lattice  Junction 


In  this  section,  the  joint  coupling  relationship  for  the  two-dimensional  rigid  joint  with 
mass  and  rotary  inertia,  shown  in  Fig.  1,  is  presented.  The  scattering  matrix  is  developed  and 
presented  both  analytically  and  graphically. 

Joint  Coupling  Matrix 

The  explicit  form  of  the  joint  coupling  matrix  for  arbitrary  two-dimensional  rigid  joints 
with  mass  and  rotary  inertia  is  given  in  [2].  From  [2],  the  joint  coupling  relationship  for  the 
joint  shown  in  Fig.  1  is 
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(31) 

0 

-men1 

0 

0 

(-1) 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

i 

0 

0 

0 

0 

(1) 

is  the  joint  coupling  i 

matrix  of  the  rigid  joint  with 

mass,  and 

H  =  {“2 

y  2 

y2 

"M2 

n2} 

(32) 

zr={«3 

>3 

^3 

(33) 

are  the  Fourier  transformed  state  vectors  at  points  2  and  3.  The  two  elements  of  the  joint 
coupling  matrix  B  which  are  enclosed  in  parentheses  are  the  modified  terms  which  ensure 


that  the  sign  conventions  for  the  joint  coupling  matrix  and  the  transfer  matrix  are  consistent. 
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The  joint  coupling  relationship  given  by  eqn.  (30)  contains  the  force  dynamic  and 
geometric  compatibility  requirements  at  the  joint.  Substituting  eqns.  (31)  through  (33)  into 


eqn.  (30)  and  writing  out  each  scalar  equation  give 

y2+I<3  =  0  (34) 

u2-y,  =  0  (35) 

¥2  "¥3  =  0  (36) 

M2-M3  =  4  w2¥2  (37) 

~V2  +N3=m(D2y2  (38) 

N2  +  V3=m(H2u2  (39) 


Eqns.  (34)  through  (36)  represent  the  compatibility  requirements  at  the  joint  and  eqns. 
(37)  through  (39)  represent  the  force  dynamic  requirements  at  the  joint.  Fig.  3  illustrates  the 
internal  forces  of  the  attached  members  acting  on  the  joint  in  accordance  with  eqns.  (34) 
through  (39). 

Scattering  Matrix 

The  process  of  transforming  eqn.  (30)  into  a  scattering  relationship  of  the  form  of  eqn. 

(1 1)  is  presented  below. 

The  state  vectors  at  points  2  and  3  in  Fig.  1  are  related  to  the  wave-mode  coordinates  at 
points  2  and  3,  respectively,  by 

(40) 

(41) 

^}2  (42) 

(43) 


Zz  =  T(C0)W2 

z  3  =  T(Q)w3 

where 

=  {w|  w+2  VV3  vvj  w~ 

W_l  =  (vv|  W+2  w;  w'3  w- 
a re  the  wave-mode  vectors  at  points  2  and  3,  re  'oectively. 

Substituting  eqns.  (40)  and  (41)  into  eqn.  (30)  gives 
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JL(fo)  0  ' 

H  0  I(CD) 


(44) 


O' 

0 


0 

0 


Substituting  eqns.  (14)  and  (31)  into  eqn.  (44)  and  multiplying  out  give 


'  0 

-1 

-1 

-1 

-1 

0 

1 

0 

0 

0 

0 

1 

0 

y32 

-4i 

-42 

0 

0 

^2-4^32 

^43  _ 4®  43 

4t3  +  4®  43 

42+4®242 

0 

0 

-(Y52  +  m(£>2) 

-(y53  +  moj2) 

Ya-m  co2 

Y52 

0 

Y6l-nm2 

0 

0 

0 

0 

-(r61  +  mco2) 

-1 

0 

0 

0 

0 

-r 

0 

0 

-1 

-l 

-l 

-l 

0 

0 

0 

"42 

-43 

43 

42 

0 

JWl 

0 

(45) 

0 

-42 

-43 

-43 

-42 

0 

kr 

0 

► 

4i 

0 

0 

0 

0 

-4, 

0 

0 

42 

43 

-43 

-42 

0 

j 

0 

From  Fig.  4,  two  wave-mode  vectors  which  represent  waves  leaving  and  entering  the 
joint  may  be  written  as 


— T 
WjOUT 

=  (w2 1 

wT2 

vv53 

(46) 

— T 
W.IN  " 

wr2 

vv23 

w33 

wT2 

w3;> 

(47) 

A  wave-mode  coordinate  wnf  represents  a  type  i  wave  (i=l,2,3)  which,  given  a  plus  (+) 
sign,  travels  along  increasing  x  or,  given  a  minus  (-)  sign,  travels  along  decreasing  x ,  at 
point  n.  For  the  system  shown  in  Fig.  1,  n  may  equal  1  through  4.  For  example,  the 
wave-mode  coordinate  represents  a  type  two  wave  (which  is  a  mode  1  flexural  wave)  at 
point  3  which  is  propagating  toward  point  4. 

Rearranging  the  order  of  the  scalar  equations  of  eqn.  (45)  gives 
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0 

-1 

-1 

0 

0 

-1 

1 

0 

0 

-1 

-1 

0 

0 

—  y32 

-^33 

-r» 

0 

0 

^2  +  4“2y32  ^3 

+/„o)2y33 

-*43 

-^42 

0 

0 

Y52-m(£>2  Y 

53 -mco2 

0 

0 

Y» 

(y61+/wco2) 

0 

0 

^53 

^52 

0 

0 

-1 

-1 

0 

0 

1 

0 

0 

-1 

-1 

0 

y* 

^33 

y32 

0 

Ya2-I„(d%2 

~Y43 

-y42 

0 

-(Y52  +  m(£>2) 

0 

0 

Ya-m  a 

i2  o 

0 

~Y» 

—  ^52 

i.  (48)  is  equivalent  t 


vv2j 

wT2 

w2 3 
^OW]  -—x+ 

w33 


- Jn+'I 

i  w2i 

d  n 

r+  [  —IN] — x-  ’ 

>3  Wh 

\*2  wT2 

1^57 


where 


0 

-1 

-1 

0 

0 

1 

0 

0 

-1 

-1 

0 

-^32 

-^33 

-y33 

-^32 

0 

ya2+L<»%2 

y43+/„co2y33 

-^43 

-^42 

0 

y52-m  co 2 

y53-m  co2 

0 

0 

(y61+ma)2) 

0 

0 

^53 

^ 52 

17 


'  0 

-1 

-1 

0 

0 

-n 

1 

0 

0 

-1 

-1 

0 

0 

Y32 

^33 

^33 

T32 

0 

Hjs  ~ 

0 

Y«-L 

~Y« 

-Y*2 

0 

0 

-(T52  +  mco2) 

-(YS3  +  mc o2) 

0 

0 

-Y6 1 

y61-mo)2 

0 

0 

-^53 

-^52 

0 

J 

Rearranging  eqn.  (49)  into  the  form  of  eqn.  (11)  gives 


w2l 

wT2 

►  =  (CD)  - 

vv3  3 

wT2 

*  J 

vv37 

1 J 

(51) 


(52) 


where 


S  (0))  -  -Djam  Hm 


(53) 


is  the  scattering  matrix  of  the  rigid  joint  with  mass  shown  in  Fig.  1.  Eqn.  (52)  is  the 
description  of  the  scattering  at  the  joint  shown  in  Fig.  1  in  terms  of  the  wave-mode 
coordinates.  The  matrix  D^m  and  the  scattering  matrix  5(co)  are  determined  symbolically 
using  MACSYMA.  The  scattering  matrix  has  the  form 


rSn 

5*12 

$13 

Su 

$ 15 

C  A 
°16 

S21 

S22 

$ 23 

$24 

$25 

$26 

5(to)  = 

S31 

532 

533 

$34 

$35 

$36 

Sn 

$42 

$43 

$44 

$45 

$46 

S5l 

$52 

$ S3 

$54 

$55 

$56 

Is- 

$62 

$63 

$64 

$65 

$66y 

The  upper-left  and  lower-right  3x3  matrices  of  £(co)  given  by  eqn.  (54)  contain  the 
reflection  coefficients  for  the  waves  which  strike  the  rigid  joint  from  point  1  and  from  point 
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4,  respectively.  The  upper-right  and  lower-left  3x3  matrices  of  5(c o)  given  by  eqn.  (54) 
contain  the  transmission  coefficients  for  the  waves  which  strike  the  rigid  joint  from  point  4 
and  from  point  1,  respectively. 

Since  member  2  is  semi-infinite  and  free  of  externally  applied  loads,  there  are  no  waves 
entering  the  joint,  at  point  3,  from  point  4.  Therefore,  substituting  eqn.  (54)  into  eqn.  (52)  and 
rewriting  give 


The  upper  3x3  matrix  of  £(ou)  in  eqn.  (55)  contains  the  reflection  coefficients  which  the 
incident  waves  at  point  2  encounter  as  they  strike  the  joint.  The  lower  3x3  matrix  of  S(c o)  in 
eqn.  (55)  contains  the  transmission  coefficients  which  the  incident  waves  at  point  2  encounter 
as  they  strike  the  joint.  From  eqn.  (55),  a  longitudinal  wave  may  reflect  and/or  transmit  a 
longitudinal  and/or  flexural  wave,  and  a  flexural  wave  may  reflect  and/or  transmit  a 
longitudinal  and/or  flexural  wave.  Systems,  such  as  the  "L"  lattice  junction,  which  exhibit 
wave-mode  conversion  between  longitudinal  and  flexural  waves  display  a  phenomenon  we 
choose  to  call  "wave-mode  conversion  of  the  second  kind".  It  is  manifested  through  the  terms 
S12,  S13,  Sa,  S2I,  S3l,  S32,  541,  Sw  5SI,  5,3,  SM,  and  Sa  of  the  scattering  matrix  S(co).  Systems  which 
exhibit  wave-mode  conversion  only  between  mode  1  and  mode  2  flexural  waves  display 
"wave-mode  conversion  of  the  first  kind”. 

From  MACSYMA  [6],  the  entries  of  the  scattering  matrix  given  in  eqn.  (55)  are 

SN 

Sij=-T  0'  =  1»2 . 6,/ =  1,2,3)  (56) 

•jo 

where 


19 


S?l  =  -[(Yn-Yjlum2v>6+ 

2[Y32Y53(Y33  -  Y32)  -  Y33Y52(Y33  -  Y32)y„m(£,A  + 

2[y43(y33  ~  Y32)  -  YA2(Y33  -  F32)]m  V  + 

[44  -  2y32y33y52y53  -  y2  (y33  -  yJ + y^yj&y + 

2[y32y43y53 + yA2y53(y33  -  2  y32) + y43yS2(y32  -  iy33) + y33y42y52]mco2 + 

2[K33y52(y43y52  -  y42y53) + y32y53(  yA2yS3  -  y43y52) + 

Ylwjj* 3  -  y32)  -  y43(y33  -  y32))  -  y33K42y52y53]]  (57) 

Sn = [2[y32y53(y43  -  y42)  -  y32yS2(y43  -  y42)]mco2+ 

2[y42y52(y33yJ2  -  y53(y33 + y32)) + y32y53y61(y43  -  y42)  - 
r52y6i(y43  -  r42) + iyy4n  (58) 

4  =  [2(y33y53(y43  -  y42)  -  y33y52(y43  -  y42))^  co2 + 

2[y43y52(y33y52  -  y53(y33 + y32)) + y33y53r61(y43  -  y42)  - 

W..^  -  ^  +  Y32YA3Yl3]]  (59) 

4  =  2[y61(y32y43  -  y33y42W + 

tfoy*  -  ym  -  y53y61(y32y43  -  y33y42)]  m 

S22  =  -[(Xi3-Y^I„m2(06+ 

2[y32(y32y53  -  y33y52) + y61(y23  -  y22)]/„ma)4 + 2y33(y43  -  y4>  V + 

[4(4  -  4) + 2y32y61(y32y53  -  y33y52)  -  y22y2]/„G)2 + 

2[2y33y61(y43  -  y42) + y33y42y53  -  y32y43y52]mco2 + 

ttOfl  -  ^42) + rn(TJjr*  -  y*yA3y52) + y43y52(y32y53  -  y33y52)]]  (6i) 

4 = 42y33(y33  -  y32)/„m  V + 

[4y33y61(y33  -  y32) + 2y33(y32y53  -  y33y52)]/„mco4 + 2y33(y43  -  y42)m  V + 

2[y33y6i(y33  —  y32) + ^33^6i(^32^53 — 

2[2y33y61(y43  -  y42) + yA3y53{y33  -  y32) + y33y42y53  -  y33y43y5>co2 + 

2[K43y53(y32y53  -  y33y52) + y61(yA3y53(y33  -  y32) + y33y42y53  -  y33y43y52) + 

y33y2(y43-y42»]  (62) 
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(63) 


S3l  =  -2[Y6l(Y32Y43  -  Y33Y42)m(>?  + 

^52^6l(^33^42  —  ^32^«)  +  ^61  (^32^43  ~  ^33^4^)] 


S?2  =  [2Y32(Y33-Y32)Izzm2U>6  + 

2[2y32y61(y33  -  y32) + y32y53  -  y^mi^ + 2  y32(y43  -  y42)m2  c»4 + 
2[y32y62(y33  -  r*) + v«(  -  v52)]4®2 + 

2[2YnY6l(Y43  —  Y42)  +  y 32^42^53  ^42^52(^33  ~  ^32)  —  ^32^43^5 + 

2[y42y52(y32y53  -  y33y52) + y42y52y61(y33  -  y32) + y32y61(y42y53  -  y45y52) + 


Y,2Y261(Y43-Ya2)]]  (64) 

Sl  =  [{Yl3-Y232)Ittm2a6  + 

2[y61(y323  -  n22) + y33(y32y53  -  y33y52)]/„mco4 + 2y32(y43  -  r>V + 

[^33^52(^33^52 — 2y32y53 — 2y33y61) + y32y53(y32^53  2y33y52) + y61(y33  -  Y32i\i„(&  + 
2[y32y43(2y61  -  y52) + y42(y33y53  -  2y32y61)]mG>2+ 

2[y42n3(y32y53  -  y33y52) + y61(y33y42yS3  -  y32y43y52) + 


Y32Yl(Y43-Y42)l]  (65) 

sHi = -[2y32y61(y33  -  Y3$„m<s? + 

2^ 61  (^ 32^53  ~  ^33^52)  +  ^32^6l(^33  —  ^32)]^k®  + 

2[y42y61(y33  -  2  y32) + y32y43y61]m(o2 + 

2[2y32y42y61(y53 — y51) + y32y43y6i(y6i —  ^52) + ^^^siO^i  ~  y52)]\  (66) 

S?2  =  [2Y32(Y43-Y42)mW+ 

2[y32y42(y53 + yJ2) + 2y32y6](y43 — y42)  -  y52(^33^42 + Y3^Y43)]mo> 

2^32^53061  ~  ^5l)  +  ^33^42^52(^52  “  ^6l)  +  ^32^6l(^43  _  ^42)H  (62) 

s43  =  [2y33(y43  -  r>V + 

2[2y33y61(y43  -  y42) + y33y52(y42  -  2  y43) + y32y43y53]mco2 + 

[2^32^43^53(^61 —  ^52) + 2r33r i3Y52(YS2 — 2yj + 2y33y61(y42y52 + 

V«-WD  (68) 


S^  =  [2Y33Y6l(Y33-Y32)Iam(£>4  + 
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(69) 


2{Y33^6l0^32^53  ^33^6l)  ^33^6l(^33^52 +  ^32^6l)]^iiW  + 

2[r43r6l(2Y33  -  y32)  -  y33y42y61]mo)2 + 

2[2y33y43y61(y61  -  y52) + y32y43y61(y53  -  y61) + WA  -  yj]\ 
5j=-[2y32(y43-y4>V+ 

2[2y32y61(y43  -  rj + y32y53(2y42  -  y43)  -  y33y42y52]mw2 + 

2[Y33Y42Y52(Y53  ~  Y6l)  +  Y32Y42Y53(2Y6l  -  Y53)  +  Y32Y6l(Y6iY43  ~ 

Y„Y53-Y42Y61)]]  (70) 

S53  =  -[2y33(y43  -  y42)mV + 

2[2y33y61(y43  -  y42) + ymy3  2  -  y33) + y33(y42y53  -  y43y52)W + 

2[y33y43y52(y53  -  y61) + y32y43yS3(y61  -  y53) + vw  v53  - 
y43y53+y43y6X-y42y61)]]  (71) 

s£ = 2[y33y42y6I(y52  -  y53) + y32yMy53  -  y52)i  (72) 

s£ = -[2[y33y52(y32  -  y33)  -  ly^s*  -  Y33)]iam  o>4 + 

2^53^1^33  -  Y32)  -  Y33Y52Y6l(Y33  -  Yn)  +  Y33Y52(Y33Ys2  -  2Y32y53)  +  Y322YS\]I„G>2  + 

2[y43y52(y32  -  2y33) + y42y52(2y33  -  y32) + y32y53(x43  -  y42)]m  co2 + 

2[y33y522(2y43  ~  y42)  ~  y32i/52^53(2y43 — y42) — y33y42y 5  2y53 + y32y53y6i(y43 — 42) — 

-  r*) + y4A52y61(2y33  -  y32) + y32y42y523]]  (73) 

5S = -[2[y32y53(y33  -  y32)  -  iw**  -  w**  w4 + 

2[y32y53y61(y33  -  14)  -  y33y52y61(y33  ~  y32) + y33y52{y33y52  ~  2  y32y43) + y32y53]4w  + 
2[y33^52(n2  -  K43) + y42y53(y33  -  2  y32)  -  y43y53(y33  -  2y32)]mco2 + 

2[Y43Y52Y53(Y33 —  ^32)  +  Y33Y52(Y43Y52  —  2Y42Y53)  -  Y33Y52Y6l(Y43  —  Y42 )  — 

Y43y53y6i(Y33  -  2  y32) + y42yS3y61(y33  -  2  y32)  -  y«yicr«  -  2  yyn  (74) 

and 

^  =  [0'33->,32)24^V  + 

spywi*  -  n3) + v53(^33  -  ^2) + 4,(^3  -  yjvarnu? + 

2[y42(y32  -  ^33) + y43(y33  -  y32))^v + 
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[2^32^53^61  (^33  ^32)  hlY SlY b\(Y V,  ^32)  2732^33^52^53  + 


Yltf2  +  Y261(Y„-Y22fya(»2  + 

2[7437S2(y32  -  2733)  +  YMY„  -  2Yn)  +  2Y„Y61(Yn  -  732)  - 


2742761(733  -  ^2)  +  ^  + 

2[iyy4  -  Vasay^ + iyy + iyy£ + 


^42^53^61(^33  2732)  +  743752^61(^32  2733)  +  Y^Y^Y^Y „  +  733742752761  + 


ViOTsS  -  *«>  -  V6.(^33  -  ^32)]] 


(75) 


The  entries  of  the  scattering  matrix  given  by  eqns.  (56)  through  (75)  are  frequency 
dependent.  Therefore,  the  system  shown  in  Fig.  1  exhibits  dispersive  scattering  as  waves  are 
reflected  at  and  transmitted  through  the  rigid  joint  with  mass  and  rotary  inertia. 

The  reflection  and  transmission  coefficients  given  by  eqns.  (56)  through  (75)  are,  in 
general,  complex  quantities.  Using  the  properties  of  a  tubular  lattice  member  and  a  rigid  joint 
with  mass  given  in  Table  1 ,  the  magnitude  and  phase  of  each  reflection  and  transmission 
coefficient  are  presented  below  as  a  series  of  plots.  These  properties  approximate  the 
hardware  currently  under  development  by  McDonnell  Douglas  and  NASA  which  will  be 
used  in  the  space  station  scheduled  for  orbital  construction  in  the  mid  1990’s. 

The  magnitude  is  calculated  as  the  square  root  of  the  sum  of  the  squares  of  the  real  and 

t" 

imaginary  parts.  The  phase  is  calculated  as  the  arctangent  of  the  imaginary  part  divided  by 
the  real  part.  The  arctangent  function  used  in  this  calculation  gives  a  value  between  ±180 
degrees.  As  a  consequence,  the  phase  plot  may  be  read  like  the  surface  of  a  cylinder  with  its 
axis  parallel  to  the  abscissa,  and  the  180  degree  and  the  -180  degree  points  on  the  ordinate 
coincident.  The  phase  of  a  reflection  coefficient  represents  a  phase  shift  between  the  reflected 
wave  and  the  incident  wave.  If  positive,  it  represents  a  phase  lead  of  the  reflected  wave  over 
the  incident  wave;  if  negative,  a  phase  lag.  Likewise,  the  phase  of  a  transmission  coefficient 
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represents  a  phase  shift  between  the  transmitted  wave  and  the  incident  wave.  If  positive,  it 
represents  a  phase  lead  of  the  transmitted  wave  over  the  incident  wave;  if  negative,  a  phase 
lag. 

The  physical  joint  is  a  spherical  orb  which  is  machined  to  accommodate  the  attachment 
of  the  lattice  members.  A  loss  of  approximately  one-third  of  its  solid  mass  occurs  due  to 
machining.  The  joint  mass  and  mass  moment  of  inertia  given  in  Table  1  serve  as  nominal 
values.  The  following  plots  present  three  curves  which  correspond  to  three  values  of  joint 
mass  and  mass  moment  of  inertia.  In  each  case,  the  joint  diameter  is  assumed  to  be  the  same, 
0.10  meter.  Therefore,  the  mass  and  mass  rotary  inertia  are  linearly  related  [7]. 

From  eqn.  (55),  it  should  be  noted  that  the  longitudinal  and  flexural  waves  are  naturally 
coupled  as  a  consequence  of  the  perpendicular  geometry  of  the  system.  For  the  properties  of 
a  tubular  lattice  member  given  in  Table  1,  the  cutoff  frequency  coc  [1,8]  of  the  mode  2 
flexural  wave  is  approximately  105,000  radians  per  second.  The  cutoff  frequency  coc  is  a 
frequency  below  which  the  mode  2  flexural  wave  is  an  evanescent  wave  and  above  which  the 
mode  2  flexural  wave  is  a  propagating  wave.  The  specific  value  of  the  cutoff  frequency  is 
strongly  dependent  on  the  shear  related  properties  of  the  lattice  member,  as  illustrated  by  eqn. 
(57)  in  [1].  Although  it  may  not  be  predominate,  the  effect  of  the  cutoff  frequency  is  apparent 
in  all  of  the  reflection  and  transmission  coefficient  plots.  In  addition,  for  frequencies  0  <  to  < 
coc,each  magnitude  curve  exhibits  a  peak  or  a  valley.  It  can  be  shown,  by  calculations  not 
shown  here,  that  the  magnitude  and  frequency  at  which  the  phenomenon  occurs  is  sensitive 
to  the  joint  mass  moment  of  inertia  and  the  shear  related  properties  of  the  lattice  members.  A 
physical  interpretation  of  this  phenomenon  and  the  possible  existence  of  a  characteristic 
equation  which  identifies  the  frequency  at  which  it  occurs  is  unclear  at  this  time  and  suggests 
further  study. 
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Figs.  5a  and  5b  give  the  magnitude  and  phase  of  the  scattering  matrix  reflection 
coefficient  Sn.  The  reflection  coefficient  Su  represents  the  amplitude  ratio  of  reflected 
longitudinal  waves  due  to  incident  longitudinal  waves  at  point  2  for  various  frequencies, 
reference  Fig.  1.  As  the  frequency  to  approaches  zero,  the  magnitude  of  S„  approaches  unity 
and  the  phase  approaches  zero  degrees.  This  suggests  that  a  longitudinal  wave  propagating 
along  member  1,  toward  the  joint,  will  perceive  the  junction  as  a  free-end  boundary  condition 
and  the  wave  will  be  reflected  accordingly.  As  the  frequency  (0  approaches  infinity,  the 
magnitude  of  Sn  approaches  unity  and  the  phase  approaches  -180  degrees.  Figs.  5a  and  5b 
illustrate  the  dynamics  between  incident  longitudinal  waves  and  reflected  longitudinal  waves 
at  point  2. 

Figs.  6a  and  6b  give  the  magnitude  and  phase  of  the  scattering  matrix  reflection 
coefficient  The  reflection  coefficient  Sn represents  the  amplitude  ratio  of  reflected 
longitudinal  waves  due  to  incident  mode  1  flexural  waves  at  point  2  for  various  frequencies, 
reference  Fig.  1.  As  the  frequency  co  approaches  zero,  the  magnitude  of  Sn  approaches  zero 
and  the  phase  approaches  90  degrees.  As  the  frequency  CO  approaches  infinity,  the  magnitude 
of  Sn  approaches  zero  and  the  phase  approaches  -180  degrees.  Figs.  6a  and  6b  illustrate  the 
dynamics  of  wave-mode  conversion  between  incident  mode  1  flexural  waves  and  reflected 
longitudinal  waves  at  point  2. 

Figs.  7a  and  7b  give  the  magnitude  and  phase  of  the  scattering  matrix  reflection 
coefficient  S13.  The  reflection  coefficient  S13  represents  the  amplitude  ratio  of  reflected 
longitudinal  waves  due  to  incident  mode  2  flexural  waves  at  point  2  for  various  frequencies, 
reference  Fig.  1.  As  the  frequency  0)  approaches  zero,  the  magnitude  and  phase  of  S13both 
approach  zero.  As  the  frequency  co  approaches  infinity,  the  magnitude  and  phase  of  S!3 
approach  zero.  Figs.  7a  and  7b  illustrate  the  dynamics  of  wave-mode  conversion  between 
incident  mode  2  flexural  waves  and  reflected  longitudinal  waves  at  point  2. 
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Figs.  8a  and  8b  give  the  magnitude  and  phase  of  the  scattering  matrix  reflection 
coefficient  S21.  The  reflection  coefficient  521  represents  the  amplitude  ratio  of  reflected  mode  1 
flexural  waves  due  to  incident  longitudinal  waves  at  point  2  for  various  frequencies, 
reference  Fig.  1.  As  the  frequency  0)  approaches  zero,  the  magnitude  of  S2i  approaches  0.5 
and  the  phase  approaches  90  degrees.  As  the  frequency  03  approaches  infinity,  the  magnitude 
of  S21  approaches  zero  and  the  phase  approaches  -180  degrees.  Figs.  8a  and  8b  illustrate  the 
dynamics  of  wave-mode  conversion  between  incident  longitudinal  waves  and  reflected  mode 
1  flexural  waves  at  point  2. 

Figs.  9a  and  9b  give  the  magnitude  and  phase  of  the  scattering  matrix  reflection 
coefficient  Sn.  The  reflection  coefficient  Sn represents  the  amplitude  ratio  of  reflected  mode  1 
flexural  waves  due  to  incident  mode  1  flexural  waves  at  point  2  for  various  frequencies, 
reference  Fig.  1.  As  the  frequency  G)  approaches  zero,  the  magnitude  of  approaches 
approximately  0.7071  and  the  phase  approaches  -135  degrees.  As  the  frequency  co 
approaches  infinity,  the  magnitude  of  Sa approaches  unity  and  the  phase  approaches  -180 
degrees.  Figs.  9a  and  9b  illustrate  the  dynamics  between  incident  mode  1  flexural  waves  and 
reflected  mode  1  flexural  waves  at  point  2. 

Figs.  10a  and  10b  give  the  magnitude  and  phase  of  the  scattering  matrix  reflection 
coefficient  Sn.  The  reflection  coefficient  represents  the  amplitude  ratio  of  reflected  mode  1 
flexural  waves  due  to  incident  mode  2  flexural  waves  at  point  2  for  various  frequencies, 
reference  Fig.  1.  As  the  frequency  c o  approaches  zero,  the  magnitude  of  SM  approaches 
approximately  0.7071  and  the  phase  approaches  -135  degrees.  As  the  frequency  G) 
approaches  infinity,  the  magnitude  of  5a  approaches  2.0  and  the  phase  approaches  -180 
degrees.  Figs.  10a  and  10b  illustrate  the  dynamics  of  wave-mode  conversion  between 
incident  mode  2  flexural  waves  and  reflected  mode  1  flexural  waves  at  point  2. 

Figs.  1  la  and  lib  give  the  magnitude  and  phase  of  the  scattering  matrix  reflection 
coefficient  S„.  The  reflection  coefficient  S31  represents  the  amplitude  ratio  of  reflected  mode  2 
flexural  waves  due  to  incident  longitudinal  waves  at  point  2  for  various  frequencies. 
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reference  Fig.  1.  As  the  frequency  0)  approaches  zero,  the  magnitude  of  Su  approaches  0.5 
and  the  phase  approaches  -90  degrees.  As  the  frequency  CD  approaches  infinity,  the  magnitude 
and  phase  of  Sn  approach  zero.  Figs.  1  la  and  1  lb  illustrate  the  dynamics  of  wave-mode 
conversion  between  incident  longitudinal  waves  and  reflected  mode  2  flexural  waves  at  point 
2. 

Figs.  12a  and  12b  give  the  magnitude  and  phase  of  the  scattering  matrix  reflection 
coefficient  S32.  The  reflection  coefficient  S32  represents  the  amplitude  ratio  of  reflected  mode  2 
flexural  waves  due  to  incident  mode  1  flexural  waves  at  point  2  for  various  frequencies, 
reference  Fig.  1.  As  the  frequency  co  approaches  zero,  the  magnitude  of  S32  approaches 
approximately  0.7071  and  the  phase  approaches  135  degrees.  As  the  frequency  CD  approaches 
infinity,  the  magnitude  S32 approaches  zero  and  the  phase  approaches  -180  degrees.  Figs.  12a 
and  12b  illustrate  the  dynamics  of  wave- mode  conversion  between  incident  mode  1  flexural 
waves  and  reflected  mode  2  flexural  waves  at  point  2. 

Figs.  13a  and  13b  give  the  magnitude  and  phase  of  the  scattering  matrix  reflection 
coefficient  S33.  The  reflection  coefficient  S33  represents  the  amplitude  ratio  of  reflected  mode  2 
flexural  waves  due  to  incident  mode  2  flexural  a .  .*  point  2  for  various  frequencies, 

reference  Fig.  1.  As  the  frequency  co  approaches  zero,  the  magnitude  of  ^approaches 
approximately  0.7071  and  the  phase  approaches  135  degrees.  As  the  frequency  co  approaches 
infinity,  the  magnitude  5„  approaches  unity  and  the  phase  approaches  zero.  Figs.  13a  and  13b 
illustrate  the  dynamics  between  incident  mode  2  flexural  waves  and  reflected  mode  2  flexural 
waves  at  point  2. 

Figs.  14a  and  14b  give  the  magnitude  and  phase  of  the  scattering  matrix  transmission 
coefficient  S41.  The  transmission  coefficient  541  represents  the  amplitude  ratio  of  transmitted 
mode  2  flexural  waves  at  point  3  due  to  incident  longitudinal  waves  at  point  2  for  various 
frequencies,  reference  Fig.  1.  As  the  frequency  co  approaches  zero,  the  magnitude  of  S41 
approaches  approximately  1.120  and  the  phase  approaches  approximately  -27  degrees.  As  the 
frequency  co  approaches  infinity,  the  magnitude  of  St!  approaches  zero  and  the  phase 
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approaches  -90  degrees.  Figs.  14a  and  14b  illustrate  the  dynamics  of  wave-mode  conversion 
between  incident  longitudinal  waves  at  point  2  and  transmitted  mode  2  flexural  waves  at 
point  3. 

Figs.  15a  and  15b  give  the  magnitude  and  phase  of  the  scattering  matrix  transmission 
coefficient  S42.  The  transmission  coefficient  S42 represents  the  amplitude  ratio  of  transmitted 
mode  2  flexural  waves  at  point  3  due  to  incident  mode  1  flexural  waves  at  point  2  for  various 
frequencies,  reference  Fig.  1.  As  the  frequency  to  approaches  zero,  the  magnitude  of  S42 
approaches  approximately  0.7071  and  the  phase  approaches  135  degrees.  As  the  frequency  to 
approaches  infinity,  the  magnitude  of  S42  approaches  zero  and  the  phase  approaches  90 
degrees.  Figs.  15a  and  15b  illustrate  the  dynamics  of  wave-mode  conversion  between 
incident  mode  1  flexural  waves  at  point  2  and  transmitted  mode  2  flexural  waves  at  point  3. 

Figs.  16a  and  16b  give  the  magnitude  and  phase  of  the  scattering  matrix  transmission 
coefficient  S43.  The  transmission  coefficient  S43  represents  the  amplitude  ratio  of  transmitted 
mode  2  flexural  waves  at  point  3  due  to  incident  mode  2  flexural  waves  at  point  2  for  '-arious 
frequencies,  reference  Fig.  1.  As  the  frequency  CO  approaches  zero,  the  magnitude  of  S43 
approaches  approximately  0.7071  and  the  phase  approaches  45  degrees.  As  the  frequency  to 
approaches  infinity,  the  magnitude  of  S43  approaches  zero  and  the  phase  approaches  -90 
degrees.  Figs.  16a  and  16b  illustrate  the  dynamics  between  incident  mode  2  flexural  waves  at 
point  2  and  transmitted  mode  2  flexural  waves  at  point  3. 

Figs.  17a  and  17b  give  the  magnitude  and  phase  of  the  scattering  matrix  transmission 
coefficient  Ssr  The  transmission  coefficient  5S1  represents  the  amplitude  ratio  of  transmitted 
mode  1  flexural  waves  at  point  3  due  to  incident  longitudinal  waves  at  point  2  for  various 
frequencies,  reference  Fig.  1.  As  the  frequency  co  approaches  zero,  the  magnitude  of  5S1 
approaches  approximately  1.120  and  the  phase  approaches  approximately  27  degrees.  As  the 
frequency  0)  approaches  infinity,  the  magnitude  of  S51  approaches  zero  and  the  phase 
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approaches  -90  degrees.  Figs.  17a  and  17b  illustrate  the  dynamics  of  wave-mode  conversion 
between  incident  longitudinal  waves  at  point  2  and  transmitted  mode  1  flexural  waves  at 
point  3. 

Figs.  18a  and  18b  give  the  magnitude  and  phase  of  the  scattering  matrix  transmission 
coefficient  Si2.  The  transmission  coefficient  S5Irepresents  the  amplitude  ratio  of  transmitted 
mode  1  flexural  waves  at  point  3  due  to  incident  mode  1  flexural  waves  at  point  2  for  various 
frequencies,  reference  Fig.  1.  As  the  frequency  0)  approaches  zero,  the  magnitude  of  SS2 
approaches  approximately  0.7071  and  the  phase  approaches  -45  degrees.  As  the  frequency  co 
approaches  infinity,  the  magnitude  of  Si2  approaches  zero  and  the  phase  approaches  -90 
degrees.  Figs.  18a  and  18b  illustrate  the  dynamics  between  incident  mode  1  flexural  waves  at 
point  2  and  transmitted  mode  1  flexural  waves  at  point  3. 

Figs.  19a  and  19b  give  the  magnitude  and  phase  of  the  scattering  matrix  transmission 
coefficient  5S3.  The  transmission  coefficient  SS3 represents  the  amplitude  ratio  of  transmitted 
mode  1  flexural  waves  at  point  3  due  to  incident  mode  2  flexural  waves  at  point  2  for  various 
frequencies,  reference  Fig.  1.  As  the  frequency  CD  approaches  zero,  the  magnitude  of  Ssj 
approaches  approximately  0.7071  and  the  phase  approaches  -135  degrees.  As  the  frequency 
co  approaches  infinity,  the  magnitude  of  S„  approaches  zero  and  the  phase  approaches  90 
degrees.  Figs.  19a  and  19b  illustrate  the  dynamics  of  wave-mode  conversion  between 
incident  mode  2  flexural  waves  at  point  2  and  transmitted  mode  1  flexural  waves  at  point  3. 

Figs.  20a  and  20b  give  the  magnitude  and  phase  of  the  scattering  matrix  transmission 
coefficient  S6I.  The  transmission  coefficient  S61  represents  the  amplitude  ratio  of  transmitted 
longitudinal  waves  at  point  3  due  to  incident  longitudinal  waves  at  point  2  for  various 
frequencies,  reference  Fig.  1.  As  the  frequency  cd  approaches  zero,  the  magnitude  of  S6l 
approaches  zero  and  the  phase  approaches  45  degrees.  As  the  frequency  CD  approaches 
infinity,  the  magnitude  of  Sfi  approaches  zero  and  the  phase  approaches  -90  degrees.  Figs. 

20a  and  20b  illustrate  the  dynamics  between  incident  longitudinal  waves  at  point  2  and 
transmitted  longitudinal  waves  at  point  3. 
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Figs.  21a  and  21b  give  the  magnitude  and  phase  of  the  scattering  matrix  transmission 
coefficient  S6r  The  transmission  coefficient  S62 represents  the  amplitude  ratio  of  transmitted 
longitudinal  waves  at  point  3  due  to  incident  mode  1  flexural  waves  at  point  2  for  various 
frequencies,  reference  Fig.  1.  As  the  frequency  a>  approaches  zero,  the  magnitude  of  S62 
approaches  zero  and  the  phase  approaches  approximately  -153  degrees.  As  the  frequency  go 
approaches  infinity,  the  magnitude  S62  approaches  zero  and  the  phase  approaches  90  degrees. 
Figs.  21a  and  21b  illustrate  the  dynamics  of  wave-mode  conversion  between  incident  mode  1 
flexural  waves  at  point  2  and  transmitted  longitudinal  waves  at  point  3. 

Figs.  22a  and  22b  give  the  magnitude  and  phase  of  the  scattering  matrix  transmission 
coefficient  Sa.  The  transmission  coefficient  SM  represents  the  amplitude  ratio  of  transmitted 
longitudinal  waves  at  point  3  due  to  incident  mode  2  flexural  waves  at  point  2  for  various 
frequencies,  reference  Fig.  1 .  As  the  frequency  0t)  approaches  zero,  the  magnitude  of  Sa 
approaches  zero  and  the  phase  approaches  approximately  -117  degrees.  As  the  frequency  go 
approaches  infinity,  the  magnitude  Sa  approaches  zero  and  the  phase  approaches  90  degrees. 
Figs.  22a  and  22b  illustrate  the  dynamics  of  wave-mode  conversion  between  incident  mode  2 
flexural  waves  at  point  2  and  transmitted  longitudinal  waves  at  point  3. 

At  arbitrary  frequencies,  the  dynamic  behavior  between  the  incident  and  the  reflected 
and  transmitted  waves  can  be  difficult  to  interpret  in  physical  terms.  Therefore,  studying  the 
limiting  values  of  the  reflection  and  transmission  coefficients  may  enhance  understanding 
and  confidence  in  the  analysis.  In  this  light,  from  the  series  of  plots  presented  above,  the 
following  observations  are  noted. 

•  As  the  frequency  to  approaches  zero,  the  reflection  and  transmission  coefficients  given 
by  S(to)  approach  either  zero  or  a  constant  nonzero  value.  A  zero  value  indicates  that 
reflection  or  transmission  does  not  occur  between  the  corresponding  wave  types.  A  nonzero 
value  indicates  that  static  coupling  exists  between  the  corresponding  wave  types  and  point 
locations. 


30 


•  As  the  frequency  Q)  approaches  zero,  a  longitudinal  wave  which  strikes  the  joint  along 
member  1  is  reflected  as  if  the  joint  end  of  member  1  were  a  free-end  boundary  condition.  As 
the  frequency  to  approaches  infinity,  a  longitudinal  wave  which  strikes  the  joint  along 
member  1  is  reflected  as  if  the  joint  end  of  member  1  were  a  fixed-end  boundary  condition. 

•  As  the  frequency  to  approaches  zero,  an  incident  mode  1  or  mode  2  flexural  wave,  at 
point  2,  will  not  result  in  a  reflected  or  transmitted  longitudinal  wave  into  member  1  or 
member  2,  respectively.  However,  an  incident  longitudinal  wave,  at  point  2,  will  result  in  a 
reflected  and  transmitted  mode  1  and  mode  2  flexural  wave  in  member  1  and  member  2, 
respectively. 

•  As  the  frequency  to  approaches  infinity,  each  incident  wave  will  encounter  complete 
reflection  at  the  joint  with  no  transmission.  This  behavior  is  expected  since  a  joint  with  mass 
and  rotary  inertia  should  appear  as  a  fixed-end  boundary  condition  at  ultra-high  frequencies 
regardless  of  the  number  or  orientation  of  any  adjacent  members  attached  to  the  joint. 

•  In  general,  each  wave  which  strikes  the  joint,  at  point  2,  will  result  in  the  reflection  of 
three  different  waves.  For  example,  an  incident  longitudinal  (type  one)  wave  at  point  2  will, 
in  general,  result  in  the  reflection  of  a  longitudinal  wave,  a  mode  1  flexural  (type  two)  wave 
and  a  mode  2  flexural  (type  three)  wave.  These  waves  will  then  propagate  along  member  1  in 
the  direction  of  decreasing  x,  refer  to  Fig.  1. 

•  In  general,  each  wave  which  strikes  the  joint,  at  point  2,  will  result  in  the  transmission 
of  three  different  waves,  at  point  3.  For  example,  an  incident  longitudinal  wave  at  point  2 
will,  in  general,  result  in  the  transmission  of  a  longitudinal  wave,  a  mode  1  flexural  wave  and 
a  mode  2  flexural  wave.  These  waves  will  then  propagate  along  member  2  in  the  direction  of 
increasing  x,  refer  to  Fig.  1 . 

In  summary,  over  various  frequencies,  the  dynamic  behavior  between  each  incident 
wave  and  each  reflected  and  transmitted  wave  is  diverse  and  unique.  In  this  analysis,  a 
longitudinal  (type  one)  wave,  a  mode  1  flexural  (type  two)  wave  and/or  a  mode  2  flexural 
(type  three)  wave  may  strike  the  joint  at  point  2.  In  general,  at  any  given  frequency  to,  each 


31 


type  of  wave  which  strikes  the  joint  may  result  in  three  reflected  waves  and  three  transmitted 
waves.  This  accounts  for  the  eighteen  scattering  matrix  terms  given  by  eqn.  (55)  which  are 
presented  and  illustrated  above  in  Figs.  5  through  22. 
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Frequency  Response  Functions  of 
"L"  Lattice  Junction  Assembly 

In  this  section,  the  frequency  response  of  each  state  variable  at  point  4  due  to  applied 
sinusoidal  loading  at  point  1  will  be  determined  and  presented  analytically  for  the  "L"  lattice 
assembly  shown  in  Fig.  1. 

The  wave-mode  coordinates  for  longitudinal  and  flexural  waves  generated  at  point  1  are 
given  by  eqn.  (29).  The  wave-mode  vector  at  point  1  is  related  to  the  wave-mode  vector  at 
point  2  according  to 

^2  =  ILMw !  (76) 

where  [1,4,5] 


is  the  wave-mode  propagation  matrix  for  member  1. 

Substituting  eqns.  (29)  and  (77)  into  eqn.  (76)  and  considering  only  waves  which 
propagate  toward  the  rigid  joint  with  mass  and  rotary  inertia  give 

'*?]  f  e+'GJ', 

■»2-  =e~‘\GJ?u+GJ!v)-  (78) 

■Wj  2  .e  +  v) 

Eqn.  (78)  contains  the  longitudinal  and  flexural  waves  at  point  2  which  strike  the  joint 
from  point  1 . 

Considering  only  waves  transmitted  through  the  joint,  substituting  eqn.  (78)  into  eqn. 

(55)  and  writing  out  each  scalar  equation  give,  after  some  manipulation, 

~  e  ^61^16^  N  +  e  +  GjjFy)  +  g  +  G35F v)  (79) 
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*35  =  e'J\xGi6FN  +  e~J\2(GvFu  +  G25FV)  +  e~\,(GMFM  +  GiSFv) 


*5;  =  e^S4iGi6FN  +  e~J\2{GuFM  +  GMFV)  +  e^SAi(GMFu  +  G35FV)  (81) 

^  =  0  (82) 

w5j  =  0  (83) 

^7  =  0  (84) 

Eqns.  (79)  through  (81)  are  the  wave-mode  coordinates  at  point  3  which  represent  the 

waves  generated  by  externally  applied  sinusoidal  loading  at  point  1.  The  wave-mode 
coordinates  at  point  3  are  related  to  the  wave-mode  coordinates  at  point  4  according  to 

uu  =  W(co)w3  (85) 

Substituting  eqn.  (77)  and  eqns.  (79)  through  (81)  into  eqn.  (85)  and  writing  out  each 
scalar  equation  give 

^4;  =  e~J^S6XGxfN  +  e~J^ + ^(G^n  +  G^v)  +  e^+J\,(GMFM  +  G35FV)  (86) 

*42+  =  e'J^ +  e^S52(G2AFM  +  GjFv)  +  e~^>+J^S53(G34FM  +  G35FV)  (87) 

*4;  =  e^+J%xGl6FN  +  e^+J\2(G2AFu  +  G^F  v)  +  e'^SA3(G3AFM  +  G35FV)  (88) 

w4;  =  0  (89) 

vv43  =  0  (90) 

w4>0  (91) 

The  wave-mode  coordinates  at  point  4  are  related  to  the  Fourier  transformed  state 
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variables  at  point  4  according  to 


Z4=J'(W)W4 


Substituting  eqn.  (14)  and  eqns.  (86)  through  (91)  into  eqn.  (92)  and  multiplying  out 
each  scalar  equation  give 

uA(cd)  =  e  S61G16F  n  + 


£  +g~(x’+',x,)5'  q  p . 

‘>62Kj2A  re  ‘363'j34  rM 
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(93) 


+  Fv 

y4((0)=  e~^+\ Gl6+e^+1\fi16  FN  + 

e^J%42GM  +  S53GM)  +  FM  + 

e^^Gv  +  S53G35)  +  e~^S43G35  +  Fv  (94) 

V4(co)=  +^32S51G16  +  +AV33S41G16  FN  + 

e^,Xl\Y33S42Gu  +  Y32S53Gm)  +  e'^Y33S43GM  +  e'^Y32Si2Gu  FM  + 

+  Y32S53G35)  +  e^Y33S43G3S  +  e'^Y^G^  Fv  (95) 

. 

Af4«a)=  c"'(X'*1,V4A,G,6  +  e‘(>’*'l'V4A1Gls  Fs  + 

c^*'1’W«GJ4+l'4jSSJG54)  +  c"a,y4JX41G„+e‘y2>>y42SMG24  Tu  + 
e^*^V4,S42G25+y4JS53GJ5)  +  c'!1’J'4JS4,G]S+«'>2>,y4ISnG2J  Fv  (96) 

V4(a>)  =  e'^VHs5lG,6+^VAVS!s41G,6  r„+ 

c'(*>*A)(y5jS4jGj4+ysAjCJ4)  +  c_B’l'5JS4!GM+«''a,yJ2Ss!G„  F„  + 

',h')(YsJ„Gv  +  Y^„a„) + e  ^rs A,G5!  +  e  ~'B’K3  fv  (97) 

. 
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(98) 


W4(a>)  =  e~j*'Y6lS6lG1J?N  + 

+  ?M  + 

+  ^VAV61S63G35  Fv 

Eqns.  (93)  through  (98)  represent  the  frequency  response  functions  of  each  Fourier 
transformed  state  variable  at  point  4  due  to  sinusoidally  applied  loads  at  point  1.  Given  the 
magnitude  and  frequency  of  the  applied  loads  at  point  1 ,  the  magnitude  and  phase  of  each 
state  variable  at  point  4  of  the  "L"  lattice  assembly,  shown  in  Fig.  1,  may  be  determined. 

It  is  important  to  note  that,  using  eqn.  (85),  the  wave-mode  vector  at  any  location  along 
member  2  may  be  determined  by  substituting  in  the  appropriate  wave-mode  propagation 
matrix  W((0).  Therefore,  once  the  wave-mode  vector  at  point  3  is  known,  the  frequency 
response  function  of  each  Fourier  transformed  state  variable  at  any  location  along  member  2 
may  be  determined  using  the  above  procedures. 

From  eqns.  (93)  through  (98),  the  time  domain  response  of  each  state  variable  at  point  4 
may  be  obtained  using  one  of  two  methods.  One  method  is  to  perform  contour  integration  in 
the  complex  plane  to  determine  the  inverse  Fourier  transform  of  the  state  variable.  An 
alternate  approach  is  to  utilize  the  concepts  of  digital  signal  processing  to  determine  the 
inverse  discrete  Fourier  transform  of  the  state  variable.  The  first  method  is  particularly 
favored  when  the  equations  are  of  known  tabulated  form.  For  cases  when  the  time  domain 
applied  loading  is  of  arbitrary  shape,  or  whose  transform  is  not  tabulated,  the  digital  approach 
is  recommended. 
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Conclusions  and  Recommendations 

The  objective  of  this  report  was  to  illustrate  the  concepts  and  features  of  wave-mode 
coordinate  analysis  by  studying  the  behavior  of  waves  propagating  through  an  "L"  lattice 
junction.  A  system  of  two  identical  semi-infinite  perpendicular  lattice  members  joined  by  a 
rigid  joint  with  mass  and  rotary  inertia  was  considered.  The  system  was  subjected  to 
sinusoidal  loading  at  a  point  along  one  of  the  lattice  members.  Each  lattice  member  was 
modeled  as  a  combined  longitudinal  rod  and  Timoshenko  beam. 

An  applied  sinusoidal  axial  force  TN,  shear  force  Tv  and  bending  moment  TM  along  a 
continuous  lattice  member  may  generate  up  to  three  independent  waves.  These  are:  one,  a 
longitudinal  (type  one)  wave;  two,  a  mode  1  flexural  (type  two)  wave  and;  three,  a  mode  2 
flexural  (type  three)  wave. 

As  the  frequency  CD  approaches  zero,  a  longitudinal  wave  which  strikes  the  joint  along 
member  1  is  reflected  as  if  the  joint  end  of  member  1  were  a  free-end  boundary  condition.  As 
the  frequency  CD  approaches  infinity,  a  longitudinal  wave  which  strikes  the  joint  along 
member  1  is  reflected  as  if  the  joint  end  of  member  1  were  a  fixed-end  boundary  condition. 

As  the  frequency  co  approaches  infinity,  each  incident  wave  will  encounter  complete 
reflection  at  the  joint  with  no  transmission.  This  behavior  is  expected  since  a  joint  with  mass 
and  rotary  inertia  should  appear  as  a  fixed-end  boundary  condition  at  ultra-high  frequencies 
regardless  of  the  number  or  orientation  of  any  adjacent  members  attached  to  the  joint. 

In  general,  each  wave  which  strikes  the  joint,  at  point  2,  will  result  in  the  reflection  of 
three  waves,  at  point  2.  For  example,  an  incident  longitudinal  wave  will,  in  general,  result  in 
the  reflection  of  a  longitudinal  wave,  a  mode  1  flexural  wave  and  a  mode  2  flexural  wave,  all 
of  which  will  propagate  along  member  1  in  the  direction  of  decreasing  x,  refer  to  Fig.  1.  In 
general,  each  wave  which  strikes  the  joint,  at  point  2,  will  result  in  the  transmission  of  three 
waves,  at  point  3.  For  example,  an  incident  longitudinal  wave  will,  in  general,  result  in  the 
transmission  of  a  longitudinal  wave,  a  mode  1  flexural  wave  and  a  mode  2  flexural  wave,  all 
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of  which  will  propagate  along  member  2  in  the  direction  of  increasing  x,  refer  to  Fig.  1.  In 
summary,  at  any  given  frequency  (0,  each  type  of  wave  which  strikes  the  joint  may  result  in 
three  reflected  waves  and  three  transmitted  waves  In  this  analysis,  a  longitudinal  wave,  a 
mode  1  flexural  wave  and/or  a  mode  2  flexural  wave  may  strike  the  joint  at  point  2.  This 
accounts  for  the  eighteen  scattering  matrix  terms  given  by  eqn.  (55)  which  were  presented 
and  illustrated. 

The  following  recommendations  are  presented  as  suggestions  for  future  research  : 

1  -  For  frequencies  0  <  Q)  <  G)c,  each  reflection  and  transmission  coefficient  exhibits  a 
peak  or  a  valley  in  its  magnitude  plot.  A  physical  interpretation  of  this  phenomenon  and  the 
possible  existence  of  a  characteristic  equation  which  identifies  the  frequency  at  which  each 
peak  and  valley  will  occur  is  unclear  and  suggests  further  investigation. 

2  -  As  the  frequency  to  approa_..es  zero,  static  coupling  exists  between  member  1  and 
member  2.  This  is  illustrated  by  the  nonzero  reflection  and  transmission  coefficients  near 
zero  frequency.  It  is  desirable  to  have  a  physical  interpretation  of  these  finite  nonzero  values 
of  the  reflection  and  transmission  coefficients.  Therefore,  it  is  recommended  that  further 
work  be  conducted  in  this  area  in  order  to  attain  a  better  understanding  of  these  quasi-static 
effects. 

3  -  The  wave  propagation  analysis  of  the  "L"  lattice  junction  is  an  important  step  in 
understanding  the  behavior  of  propagating  waves  in  large  space  structures.  However,  there 
are  joint  geometries  of  greater  complexity  that  must  eventually  be  addressed.  It  is 
recommended  that  a  wave  propagation  analysis  be  performed,  using  wave-mode  coordinates, 
of  a  "T"  lattice  junction. 

4  -  It  is  desirable  to  obtain  the  time  domain  response  of  each  state  variable  at  an  arbitrary 
location  along  the  "L"  lattice  assembly.  It  is  recommended  that  a  procedure  which  utilizes  the 
concepts  of  digital  signal  processing  be  pursued  in  order  to  provide  time  domain  responses 
from  the  frequency  response  functions  obtained  using  wave-mode  coordinates. 
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5  -  Finally,  it  is  recommended  that  the  results  presented  in  this  report  be  verified 
experimentally. 

Wave-mode  coordinate  analysis  is  an  effective  tool  in  studying  the  behavior  of 
propagating  waves  through  an  "L"  lattice  junction.  It  provides  a  well  structured  and  relatively 
straightforward  procedure  to  determine  the  reflection  and  transmission  coefficients  for 
longitudinal  and  flexural  waves  which  may  strike  the  joint.  In  addition,  the  frequency 
response  of  each  transformed  state  variable  may  be  determined  at  any  location  in  the  "L" 
lattice  assembly.  Unlike  finite  element  methods,  the  wave-mode  coordinate  approach  to  wave 
propagation  analysis  does  not  require  remodelling  in  order  to  increase  analysis  accuracy  or  to 
obtain  solutions  at  arbitrary  locations.  This  is  a  very  attractive  feature  which  emphasizes  the 
analysis  capabilities  and  versatility  of  wave-mode  coordinate  analysis. 
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Table  1  Geometric  and  material  properties  of  aluminum  lattice  "L"  joint  assembly. 


|  Aluminum  Alloy  Lattice  "L"  Joint  Assembly 

Tubular  Lattice  Member 

Outside  Diameter  : 

0.0750  m 

Wall  Thickness  : 

0.0125  m 

Cross-Sectional  Area  (A) : 

2.45E-03  m2t 

Section  Moment  of  Ineria  (/)  : 

1 .25E-06  m4 

Shear  Correction  Coefficient  (k)  : 

0.568 

Mass  Density  (p) : 

2.70E+03  Kg/m3 

Section  Rotary  Inertia  ( J) : 

3.37E-03  Kg-m 

Tensile  Modulus  ( E ) : 

7.00E+10  N/m2 

Shear  Modulus  (G) : 

2.60E+10  N/m2 

Rigid  Joint 

Joint  Diameter : 

0.100  m 

Mass  (m) : 

1.000  Kg  $ 

Mass  Moment  of  Inertia  (IB) : 

1.00E-03  Kg-m2* 

+  E-03  s  10 3 

$  Assumes  joint  is  spherical  and  approximately  two-thirds  solid. 
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Fig.  2a  Simulated  rigid  massless  joint  with  no  spatial  extent  at  point  of  applied 
loading  along  lattice  member. 


Fig.  2b  Wave-mode  coordinates  entering  and  leaving  simulated  rigid  massless 
joint  at  point  of  applied  loading. 
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Fig.  3  Components  of  state  vectors  at  "L"  lattice  joint  with  mass  and  rotary 
inertia. 
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Fig.  6a  Magnitude  of  scattering  matrix  reflection  coefficient  S12. 


Fig.  6b  Phase  of  scattering  matrix  reflection  coefficient  S12. 


A:  m-0.5  kg.  lzz-,0005  kg-mA2 
B:  m-1 .0  kg,  lzz-,0010  kg-mA2 
C:  m-2.0  kg.  lzz-,0020  kg-mA2 


Fig.  7a  Magnitude  of  scattering  matrix  reflection  coefficient  S]3. 
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Fig.  8a  Magnitude  of  scattering  matrix  reflection  coefficient  S21. 


Fig.  8b  Phase  of  scattering  matrix  reflection  coefficient  Sa. 
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Fig.  12a  Magnitude  of  scattering  matrix  reflection  coefficient  S32. 


Fig.  12b  Phase  of  scattering  matrix  reflection  coefficient  S 


Phase  of  S  ^  Degrees  Magnitude  of  S 


Fig.  13a  Magnitude  of  scattering  matrix  reflection  coefficient  SM. 


Fig.  13b  Phase  of  scattering  matrix  reflection  coefficient  SM 
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A:  m-0.5  kg,  lzz-.0005  kg-mA2 
B:  m-1 .0  kg,  lzz-,0010  kg-mA2 
C:  m-2.0  kg.  lzz-.0020  kg-mA2 
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Fig.  15a  Magnitude  of  scattering  matrix  transmission  coefficient  S, 


A:  m-0.5  kg,  lzz-.0005  kg-mA2 
B:  m-1.0  kg,  Izz-.OOIO  kg-mA2 
C:  m=2.0  kg,  lzz-.0020  kg-mA2 


-i - 1 - 1 - 1 - 1 - i - i - 1 - 1 - 1 - 1 - 1 - 1  i - 1  i  i  r 

20  40  60  80  100  120  140  160  180 

Frequency  -  Radians  /  Second 


Fig.  15b  Phase  of  scattering  matrix  transmission  coefficient  S, 
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Fig.  16a  Magnitude  of  scattering  matrix  transmission  coefficient  S43. 


Fig.  16b  Phase  of  scattering  matrix  transmission  coefficient  S. 


A:  m-0.5  kg,  lzz-,0005  kg-mA2 
B:  m=1 .0  kg,  Izz-.OOI 0  kg-mA2 
C:  m-2.0  kg,  lzz-,0020  kg-mA2 
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Fig.  17a  Magnitude  of  scattering  matrix  transmission  coefficient  S, 


A:  m=0.5  kg,  lzz=.0005  kg-mA2 
B:  m=1.0  kg,  lzz=.0010  kg-mA2 
C:  m=2.0  kg,  lzz=.0020  kg-mA2 
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Fig.  17b  Phase  of  scattering  matrix  transmission  coefficient  S 
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A:  m*=0.5  kg,  lzz-,0005  kg-mA2 
B:  m-1.0  kg,  lzz«.0010  kg-mA2 
C:  m=2.0  kg,  tzz-.0020  kg-mA2 


Fig.  20a  Magnitude  of  scattering  matrix  transmission  coefficient  S6r 
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A:  m=0.5  kg,  lzz-.0005  kg-mA2 
B:  m-1.0  kg,  Izz-.OOIO  kg-mA2 
C:  m-2.0  kg,  lzz-.0020  kg-mA2 


Appendix  A : 

Summary  of  Wave-Mode  Coordinate  Analysis 
of  One-Dimensional  Continuous  Lattice  Member 
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In  this  appendix,  a  summary  of  the  concepts  and  nomenclature  of  wave-mode  coordinate 
analysis  used  to  study  wave  propagation  in  continuous  one-dimensional  lattice  members  is 
presented.  The  lattice  member  is  modeled  as  an  elastic  longitudinal  rod  for  axial  motion  and 
a  Timoshenko  beam  for  flexural  motion.  A  more  detailed  discussion  of  the  major 
developments  that  follow  are  discussed  in  [Alj. 

State  Vector 

A  section  of  the  lattice  member  is  shown  in  Fig.  Al,  which  illustrates  the  local 
coordinate  system  defining  the  positive  orientation  of  axial  displacement  u,  transverse 
displacement  y,  rotation  y,  axial  force  N,  shear  force  V  and  moment  M.  The  member  is 
assumed  to  be  uniform,  with  cross-sectional  area  A,  section  moment  of  inertia  /,  tensile 
modulus  E,  shear  modulus  G,  mass  density  p  and  a  shear  correction  coefficient  K. 

The  Fourier  transformed  state  vector  at  point  1  and  point  2  are 


r(x1,CD)  =  {M1  y, 

Vi 

M,  V \ 

«,} 

(Al) 

f(x2,ta)  =  {u2  y2 

V2 

W2  V2 

(A  2) 

where  the  superscript T  indicates  a  matrix  transpose  operation.  An  overbar  signifies  a  Fourier 
transform  and  an  underbar  signifies  a  column  vector  or  matrix. 

Transfer  Matrix 

The  Fourier  transformed  state  vector  at  point  1  and  point  2  are  related  by  the  transfer 
matrix  relationship 

l(x2,(a)  =  7  (cu)  z(x,  ,co)  (A3) 

where  7(a))  is  the  transfer  matrix  of  order  n,  where  n  is  the  (even)  number  of  components  of 
the  state  vectors.  The  transfer  matrix  for  the  combined  longitudinal  rod  and  Timoshenko 
beam  lattice  member  shown  in  Fig.  Al  is 
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c0  =  A  cosh  X,  +  >4  cos  Xj) 

(A  15) 

(A*  ,  %  ,  ) 

c,  =  A  —  sinh  Ag  +  —  sin  Ag 
^^3  A2  ^ 

(A  16) 

c2  =  A(cosh  Ag  -  cos  Aj) 

(A  17) 

(  sinh  A,  sinA,^ 

C3_A(  A3  A*  J 

(A  18) 

The  eigenvalues  and  eigenvectors  of  the  transfer  matrix  T(us)  form  the  wave-mode 
propagation  matrix  and  the  wave-mode  matrix,  respectively,  of  the  lattice  member. 
Wave-Mode  Propagation  Matrix 

There  are  n  eigenvalues  of  the  transfer  matrix  T(co).  For  uniform  lattice  members,  the 
eigenvalues  occur  in  pairs,  where  the  eigenvalues  of  an  eigenpair  are  algebraic  inverses  of 
each  other.  Each  eigenpair  represents  an  independent  wave  which  may  propagate  along  the 
lattice  member.  Therefore,  there  are  n/2  possible  independent  propagating  waves  which  may 
traverse  the  lattice  member.  The  physical  interpretation  of  the  fact  that  the  eigenvalues  of 
T(wi)  occur  in  pairs  is  that  identical  waves  may  propagate  in  the  direction  of  increasing  x  or  in 
the  direction  of  decreasing  x.  The  wave-mode  propagation  matrix  is  an  n  x  n  diagonal  matrix 
whose  diagonal  elements  are  the  eigenvalues  of  Z(co).  The  wave-mode  propagation  matrix  of 
the  lattice  member  shown  in  Fig.  A1  is 
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(A  19) 


From  eqn.  (A  19),  there  are  three  eigenpairs  and,  therefore,  three  independent  propagating 
waves  which  may  simultaneously  exist  in  the  lattice  member.  The  eigenvalues  located  at 
{row,  column )  locations  (1,1)  and  (6,6)  of  W(w)  represent  longitudinal  (type  one)  waves 
which  may  propagate  along  the  lattice  member  toward  increasing  and  decreasing  x, 
respectively.  The  eigenvalues  located  at  locations  (2,2)  and  (5,5)  of  W(tt>)  represent  mode  1 
flexural  (type  two)  waves  which  may  propagate  along  the  lattice  member  toward  increasing 
and  decreasing  x,  respectively.  The  eigenvalues  located  at  locations  (3,3)  and  (4,4)  of  W(co) 
represent  mode  2  flexural  (type  three)  waves  which  may  propagate  along  the  lattice  member 
toward  increasing  and  decreasing  x,  respectively. 

Wave-Mode  Matrix 

For  each  eigenvalue  of  given  by  VF(«),  there  exists  an  eigenvector.  The 
eigenvectors  of  T(wi)  are  the  columns  of  the  wave-mode  matrix.  The  wave-mode  matrix  for 
the  lattice  member  shown  in  Fig.  A1  is 
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The  wave-mode  matrix  T(co)  exhibits  features  analogous  to  the  modal  matrix  of  discrete 
vibratory  systems  [A2J.  Each  column  in  y(co)  represents  the  normalized  modal  configuration 
of  each  Fourier  transformed  state  variable  for  each  type  of  wave. 

Wave-Mode  Vector 

A  new  vector  ~w(x, co)  is  defined  according  to 

vv(x  ,co)  =  rl  (c o)l (x  ,co)  0421) 

where  y£(x,w),  a  column  vector  of  length  n,  is  called  the  wave-mode  vector.  The  components 

of  jJvCc.co)  are  called  wave-mode  coordinates.  The  wave-mode  vector  at  point  2  is  related  to 
the  wave-mode  vector  at  point  1  according  to 

Tv(jc2,G))  =  W_  (co)vv  (x,  ,co)  ( A  22) 

Writing  out  eqn.  (A22)  gives 
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Eqn.  (A23)  is  the  diagonalized  form  of  the  transfer  matrix  relationship  given  by  eqn.  (A3). 
The  components  w;,  ^  and  represent  longitudinal,  mode  1  flexural,  and  mode  2  flexural 
waves,  respectively,  which  propagate  along  the  lattice  member  toward  increasing  x. 

Similarly,  wj,  ^  and  wj  represent  longitudinal,  mode  1  flexural,  and  mode  2  flexural  waves, 
respectively,  which  propagate  along  the  lattice  member  toward  decreasing  x. 

Finally,  from  eqn.  (A21),  it  follows  that  the  Fourier  transformed  state  vector  at  any 
specific  location  along  the  lattice  member  may  be  determined  from  the  wave-mode  vector  at 
that  location  according  to 

iQc  ,co)  =  y(co)w(x  ,to)  (A  24) 

In  summary,  the  concepts  of  wave-mode  coordinate  analysis  are  developed  in  the 
frequency  domain  through  the  Fourier  transform.  From  wave-mode  coordinate  analysis,  the 
combined  longitudinal  rod  and  Timoshenko  beam  lattice  member,  shown  in  Fig.  Al,  will 
support  as  many  as  three  independent  propagating  waves.  The  waves  that  the  lattice  member 
will  support  include  a  longitudinal  (type  one)  wave,  a  mode  1  flexural  (type  two)  wave  and  a 
mode  2  flexural  (type  three)  wave. 
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Fig.  A1  Local  coordinate  system  and  state  variables  of  lattice  member. 
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